中 华 猕 猴 桃 在 中 国 潜在 分 布 及 其 对 气候 变化 啊 应 的 研究 


GEO EK O0. ZR COO 48 RAO? ; 30] Jf O9? 


(1) 四 川 农 业 大 学 农学 院 成 都 611130; (2). 四 川 省 农村 经 济 综合 信息 中 心 ”成 都 610072; (3). 四 川 苍 溪 猕 猴 桃 研究 所 ” 苍 溪 
628400 


摘 要 : 中 华 狂 猴 桃 为 中 国 特有 果 种 ， 由 于 其 独特 的 口感 和 较 高 的 经 济 价值 ， 近 年 来 种 植 规模 逐年 扩大 。 在 引种 过 程 中 ,由 
于 缺乏 合理 的 布局 规划 和 适 生 性 分 析 ,出现 了 品种 单一 化 、 易 感 病虫害 等 问题 。 近 年 来 四 川 、 陕 西 、 贵 州 、 重 庆 和 湖北 等 狂 
猴 桃 主 产 省 份 相继 开展 了 猕猴 桃 气候 适宜 性 区 划 的 研究 , 但 目前 的 研究 多 未 考虑 未 来 气候 变化 对 猕猴 桃 种 植 分 布 的 影响 ， 且 
伴随 着 气候 变化 的 加 剧 ， 已 有 的 研究 结果 已 不 能 完全 适应 实际 生产 的 需求 本文 运 用 生态 位 模型 软件 MaxEnt， 模 拟 和 预测 气 
候 变 化 背景 下 大 尺度 范围 中 华 猕猴 桃 适 生 区 分 布 及 其 变化 的 可 行 性 ， 以 利于 科学 地 优化 产业 结构 、 促 进 产业 发 展 。 基 于 当前 
数据 和 IPCC AR5 提出 的 3 种 气候 情景 以 及 中 华 儿 猴 桃 的 分 布 信息 ,采用 MaxEnt 生态 位 模型 和 ArcGIS 预测 了 中 华 狐 猴 桃 的 
适 生 区 及 未 来 的 变化 趋势 , 用 受 试 者 工作 特征 曲线 (receiver operating characteristic curve, ROC 曲线 ) 检 测 模 型 精度 、 刀 切 法 
(jackknife test) ff t EEDA E., ARRA, 基于 当前 和 未 来 情景 构建 的 中 华 独 猴 桃 地 理 分 布 模型 的 AUC(area under 
curve) 值 均 达 到 “ 极 好 ”的 标准 ,说 明 模 型 预测 结果 可 用 于 本 研究 。 当 前 气候 条 件 下 ， 中华 狗 猴 桃 的 高 适 生 区 主要 在 四 川 、 陕 
西 、 重 庆 、 湖 北 、 贵 州 、 浙 江 、 湖 南 、 安 徽 、 河 南 、 江 苏 和 甘肃 等 省 份 , 面积 达 1.01x105km2。 中 适 生 区 则 以 高 适 生 区 为 中 
心 向 外 扩散 ,包括 河南 、 湖 北 、 安 徽 、 江 苏 和 山东 等 地 ,面积 为 6.79x105km2?。RCP2.6 和 RCP4.5 排放 情景 下 ， 中 华 独 猴 桃 高 
适 生 区 的 分 布 、 面 积 及 中 心 点 位 置 都 有 所 不 同 ,面积 均 旦 增加 趋势 ,RCP8.5 排放 情景 下 ,高 适 生 区 面积 呈 减 少 趋势 。RCP4.5 
和 RCP8.5 排放 情景 下 ,中华 独 猴 桃 高 适 生 区 中 心 点 均 有 向 北 移动 趋势 。MaxEnt 模型 对 未 来 气候 交 化 条 件 下 中 华 独 猴 桃 适 生 
区 的 准确 模拟 与 预测 具有 潜在 应 用 价值 ， 对 该 果树 的 气候 适宜 性 区 划 具 有 重要 指导 意义 。 
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Potential distribution of Actinidia chinensis in China and its predicted response to climate 
change" 
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(1. College of Agronomy, Sichuan Agricultural University, Chengdu 611130, China; 2. Sichuan Provincial Rural Economic 
Information Center, Chengdu 610072, China; 3. Kiwifruit Institute of Cangxi County, Cangxi 628400, China) 


Abstract: Kiwifruit (Actinidia spp.,) belonged to Actinidiaceae,, is a type of perennial deciduous woody liana and an important class 
of berry fruit. With rich sugar, protein, amino acids, vitamins and especially high vitamin C content, the kiwifruit is known as “the 
king of the fruit" and has a good market prospect. A. chinensis is a species endemic in China with a fast-expanding planting area due 
to its unique subtle flavor and high economic value. Optimization of planting scale and distribution of the crop has been the major 
concern for regional planning. The objective of this study was to test and determine the possibility of using the MaxEnt (the 
maximum entropy) model to simulate and predict future large-scale distribution of A. chinensis. Based on current environmental 
factors, three future climate scenarios suggested in the IPCC fifth report and current distribution sites of A. chinensis, the MaxEnt 
model was used in combination with ArcGIS to predict the potential geographic distribution and trend of change of A. chinensis in 
China. The dominant factors were chosen using the jackknife test and the receiver operating characteristic curve (ROC) used to 
evaluate the simulation. The results showed that high value of area under curve (AUC) denoted good results which significantly 
differed from random predictions. Based on the evaluation criterion, the accuracies of the predictions of A. chinensis potential 
distribution in the current and future periods were excellent. The predicted result of the MaxEnt model was imported into 
ArcGIS10.0 for further analysis and showed that under present climatic conditions, the total suitable area was 26.9296 of the total 
land area in China. The potential distribution was highly consistent with the locations of specimen records and field surveys. The 
highly suitable areas were in Sichuan, Shaanxi, Chongqing, Hubei, Guizhou, Zhejiang, Hunan, Anhui, Henan, Jiangsu and Gansu 
Provinces. The areas of highly suitable habitat in the main producing provinces were analyzed statistically. The results showed that 
under the current conditions, the most suitable area for A. chinensis cultivation was 1.01x106 km?, accounting for 38.94% of the total 
suitable areas. The moderately suitable areas were in Henan, Hubei, Anhui and Shandong Provinces, with the area of 6.79x10? km?, 
accounting for 26.26% of the total suitable areas. Comparison of future suitable areas with current suitable areas showed that areas of 
high suitability increased under scenarios RCP2.6 and RCP4.5, but decreased under scenario RCP8.5. Under scenarios RCP4.5 and 
RCPS.5, the mean center of highly suitable area of A. chinensis moved northward. The result showed that the MaxEnt model was 


highly reliable in determining not only the range of geographic distribution of A. chinensis, but also in identifying dominant 


environmental factors driving the geographic distribution. Whereas climate was a decisive factor in species distribution, change in 
distribution pattern of species was the most direct effect of climate change. The results provided a critical reference base for A. 
chinensis plantation pattern and countermeasures to cope with climate change in China. 
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concentration pathway scenario 


生物 与 气候 之 间 相 互 作用 对 物种 的 分 布 影响 巨大 , 气候 变化 不 仅 可 以 改变 物种 所 在 的 生态 系统 ， 且 造 
成 的 改变 往往 是 不 可 逆 的 。 在 区 域 尺度 上 , 气候 是 影响 物种 分 布 的 主要 环境 因素 止 。 农 业 生 产 强 烈 依赖 气 
候 生态 条 件 , 气候 变化 将 改变 农作物 已 有 的 适宜 性 区 划 ， 影 响 其 生长 、 发 育 及 产量 ?1。 物种 分 布 模型 
(species distribution models, SDM) 是 研究 实际 气候 变化 对 物种 分 布 影响 的 重要 工具 ， 近 年 来 该 类 模型 在 国 
内 外 得 到 广泛 应 用 ,目前 常用 的 有 最 大 炉 模 型 MaxEnt (maximum entropy model), 基于 遗传 算法 的 规则 组 合 
模型 GARP(genetic algorithm for rule set produciton model)、 生 物 气候 模型 BIOCLIM(biological climatic model) 
和 气候 动态 模拟 软件 Climex (match climate and compare location model) 等 约 14 种 物种 分 布 模型 571,MaxEnt 
模型 的 模拟 精度 要 高 于 其 他 模型 ,加 之 该 软件 具有 运行 时 间 短 、 操 作 简 便 、 运 行 结果 稳定 和 所 需 样 本 量 小 
等 特点 在 业内 得 到 了 广泛 认可 虽 。 近 年 来 , 国内 学 者 利用 MaxEnt 模型 成 功 研究 了 气候 变化 对 水 稻 (Oryza 
sativa) 中 、 玉 米 (Zea may) ®, /|NZE (Triticum aestivum)t 等 粮食 作物 的 影响 ， 并 分 析 了 气候 因子 与 种 植 分 布 
之 间 的 关系 , 取得 了 较 的 模拟 效果 。 

猕猴 桃 (4ctinidia spp.), 又 称奇 异 果 , 隶属 于 猕猴 桃 科 (Actinidiaceae)， 猴 桃 属 (4ctinidia)， 为 浆果 类 落叶 
滕 本 果树 ,是 20 世纪 和 人工 驯化 栽培 野生 果树 最 有 成 就 的 四 大 果树 之 一 [ 引 。 中国 是 世界 猕猴 桃 起 源 中 心 ， 猕 
猴 桃 属 66 个 种 中 有 62 个 种 原 产 于 中 国 031。 近 年 来 ,， 随 着 我 国 农 业 产业 结构 调整 ， 猕 猴 桃 产业 发 展 迅 速 , 狂 
猴 桃 的 栽培 面积 和 产量 都 得 到 长 足 的 发 展 ,目前 我 国 猕猴 桃 年 产量 已 超过 意大利 ,成 为 世界 第 一 产 出 国 。 
rH AES (Actinidia chinensis) fl Je W&AWABe BE (Actinidia deliciosa) 是 我 国 栽培 最 广泛 的 两 个 种 。 中 华 独 猴 桃 
为 中 国 特有 果 种 ， 由 于 其 口感 独特 、 经 济 价值 高 ,近年 来 种 植 规 模 逐 年 扩大 。 在 引种 过 程 中 ， 由 于 缺乏 合理 
的 布局 规划 和 适 生 性 分 析 ， 出现 了 品种 单一 化 、 易 感 病虫害 等 问题 叶 避 。 和 针对 上 述 问题 , 近年 来 四 川 041、 
陕西 U7、 贵州 09、 重 庆 [9 和 湖北 E09 等 省 份 开 展 了 狂 猴 桃 气候 适宜 性 区 划 的 研究 ， 并 在 实际 生产 中 起 到 了 各 
极 作 用 。 但 目前 的 研究 多 未 考虑 未 来 气候 变化 对 猕猴 桃 种 植 分 布 的 影响 ， 且 伴随 着 气候 变化 的 加 剧 , 已 有 
的 研究 结果 已 不 能 完全 适应 实际 生产 的 需求 ， B FUR AASETERSERT T EIAN UE E, 以 利于 
科学 地 优化 产业 结构 、 促进 产业 发 展 。 本 研究 利用 MaxEnt 模型 ， 结 合 GIS 技术 , Butt rp 4e pet nn BET] 
潜在 分 布 范 围 , Tul. OE EG. 分析 中 华 猕猴 桃 在 未 来 3 种 气候 情景 RCP2.6、RCP4.5 和 RCP8.5) 下 分 布 范 围 
和 空间 格局 的 变化 规律 ， 则 在 探索 中 华 狂 猴 桃 合理 区 划 的 有 效 手段 ， 以 期 为 其 生产 布局 及 应 对 未 来 气候 变 
化 提供 基本 的 数据 支持 。 

1 材料 与 方法 
1.1 预测 软件 来 源 

选用 基于 最 大 炉 理论 的 MaxEnt 模型 软件 来 预测 中 华 狗 猴 桃 在 中 国 的 适 生 分 布 。 登 录 普 林 斯 顿 大 学 
Robert Schapire 计算 机 科学 研究 中 心 的 MaxEnt 主页 可 免费 下 载 该 软件 ， 当 前 最 新 版 本 为 3.3.3k hic? 

选用 ArcGIS 软件 将 MaxEnt 模型 的 运算 结果 投射 到 地 图 上 显示 ,预测 中 华 狂 猴 桃 的 实际 分 布 和 潜在 分 
布 。 

1.2 物种 分 布 记录 

运用 生态 位 模型 模拟 物种 适 生 区 时 [33J， 首 先 应 获取 目标 物种 的 实际 分 布 数据 ， 本 研究 中 四 川 境内 中 
华 猕猴 桃 的 分 布 数据 为 实地 考察 获得 ,采用 GPS 记录 经 纬度 ， 其 他 地 区 的 分 布 数据 则 通过 查询 物种 分 布 数 
据 库 和 检索 国内 外 公开 发 表 的 中 华 猕 猴 桃 相关 的 期 刊 论文 获得 。 本 研究 查询 的 数据 库 包 括 “ 国 际 农业 与 生 
物 科 学 中 心 (CABL http://www.cabi.org/)” 数 据 库 、“ 全 球 物种 多 样 性 信息 库 (GBIEF, http://www.gbif.org/” 和 “ 教 
学 标本 资源 共享 平台 ”(http://mnh.scu.edu.cn/)。 通 过 上 述 方式 共 获 得 分 布点 283 个 ,使 用 Google Earth 软件 
查询 分 布点 的 经 纬度 信息 ,按照 MaxEnt 软件 要 求 ， 去 除 重复 记录 、 模 糊 记 录 和 邻近 记录 的 分 布点 。 将 最 终 
确定 的 经 纬度 数据 使 用 Excel 处 理 , 保存 格式 为 *.CSVDP4。 

153 ”环境 变量 的 选择 

构建 物种 生态 位 模型 ,还 应 选取 合适 的 环境 变量 。 本 研究 从 worldclim 网 站 下 载 了 67 个 环境 变量 , 包 
括 19 个 具 ua cu ci EX Pi E. ta 根据 个 同 需求 计算 
演 生 而 来 的 气候 变量 ,可 反映 温度 与 降水 的 特点 及 季节 性 变化 特征 )、 月 平均 温度 、 月 最 高 稳定 、 月 最 低温 


ux 


度 和 月 降水 量 。 上 述 数 据 的 空间 分 辨 2.5arc-minutes( 约 4.5 km2) 。 当 前 情景 的 数据 在 worldclim 网 站 下 载 。 
气候 变化 专门 委员 会 UPCO) 第 5 次 评估 报告 共 发 布 的 4 种 典型 浓度 路 径 P5lGepresentative concentration 
pathways, RCP), EH RCP2.6、RCP4.5、RCP6.0 和 RCP8.5。 其 中 RCP4.5 和 RCP6.0 均 为 中 等 CO» HERES 
R, 且 前 者 优先 性 大 于 后 者 ， 因 此 本 研究 选择 RCP2.6 (最 低 CO» 排放 情景 )、RCP4.5( 中 等 CO: 排放 情景 ) 和 
RCP8.5( 最 高 CO» 排放 情 景 )3 种 排放 情景 作为 未 来 气候 数据 。 未 来 时 段 包 括 2030s(2021 一 2040 年 )、2050s 
(2041—2060 年 )、2070s(2061 一 2080 年 ) 和 2080s(2071—2090 fF), 数据 在 国际 热带 农业 中 心 (International 
Centre for Tropical Agriculture, CIAT) ii FES, 
14 ”中华 猕猴 桃 生境 模型 构建 与 结果 评价 

模型 的 使 用 : 输入 中 华 猕 猴 桃 分 布点 数据 和 环境 因子 数据 ,随机 选取 75% 的 中 华 猕猴 桃 分 布点 作为 训 
练 集 (training data) 建 立 预 测 模型 剩余 25% 的 中 华 猕猴 桃 分 布点 作为 测试 集 (test data) 验 证 模型 ， 选择 刀 切 
法 (jackknife tesb 测 定 各 变量 权重 ,选择 创建 环境 变量 响应 曲线 ， 其 余 参数 均 选 择 模型 的 默认 值 27。 

模型 输出 的 数据 为 ASCI, 利用 ArcGIS 的 ArcToolbox 的 格式 转换 工具 , 将 该 数据 转 为 Raster 格 
式 , 使 该 结果 可 在 ArcGIS 中 显示 。MaxEnt 软件 模拟 输出 的 结果 值 在 0~1 之 间 , 值 越 接 近 1 表示 物种 越 可 
能 存在 。 本 研究 采用 受 试 者 工作 特征 曲线 (ROC 曲线 , receiver operating characteristic curve) 分 析 法 对 模拟 结 
有 果 进 行 精度 评价 ,曲线 下 面积 值 (area under curve, AUC) 反 映 预 测 精度 ,理论 上 AUC 指标 取 值 范围 为 0.5~1， 
EE T. 具体 评价 标准 见 表 108291, 

表 1 曲线 下 面积 值 (AUC 值 ) 取 值 范围 及 其 与 MaxEnt 模型 准确 性 的 关系 


Table 1 Relationship between area under curve (AUC) and the accuracy of the MaxEnt model 


AUC 取 值 范围 评价 标准 
Range of AUC value Evaluation criterion 
0.5 X AUC«0.6 失败 Fail 
0.6: AUC«0.7 较 差 Poor 
0.7<AUC<0.8 一 般 Fair 
0.8 AUC«0.9 好 Good 
0.9: AUC«1.0 极 好 Excellent 


1.5 中华 猕猴 桃 适 生 等 级 划分 

MaxEnt 模型 输出 结果 为 ASC 了 [格式 文件 ,首先 使 用 ArcGIS 的 ArcToolbox 的 格式 转换 工具 ,将 该 数据 
转 为 Raster 格式 , 使 该 结果 可 在 ArcGIS 中 显示 , 而 后 使 用 “提取 分 析 ” 功 能 得 到 中 华 猕 猴 桃 在 研究 区 的 存在 
概率 分 布 图 2430。MaxEnt 软件 模拟 输出 的 结果 值 为 0~1, 值 越 接 近 1 表示 物种 越 可 能 存在 。 参 考 IPCC 报 
告 B139 关 于 评估 可 能 性 的 划分 方法 ,结合 中 华 狂 猴 桃 的 实际 情况 ， 利 用 “Reclassify” 功 能 ， 划 分 分 布 值 等 级 
及 相应 分 布 范围 ， 并 使 用 不 同 颜 色 表 示 , 划分 标准 为 : 存在 概率 <0.05 为 不 适 生 区 ; 0.05 三 存在 概率 <0. 33 为 
低 适 生 区 ; 0.33 科 存在 概率 <0.66 为 中 适 生 区 ; 存在 概率 之 0.66 为 高 适 生 区 。 
L6 ”几何 中 心 及 位 移 计 算 

参考 Yue 等 B31 的 计算 方法 , 统计 了 高 适 生 区 在 未 来 不 同时 期 的 面积 变化 情况 以 及 其 中 心 点 的 位 移 , 具 
体 计算 公式 如 下 : 


x(t) = y SEU) 


il S(t) 
I i e Yi 
yA) = 2, e NE c e 


(4) 

式 中 :7 表示 不 同 的 时 段 ; 7 为 高 适 生 区 的 单位 栅 格 数目 ; si:(D) 为 1 时 段 第 ii=0, 1, 2, .…, DD 个 单位 栅 格 的 面积 ; 
SG 为 ! 时 段 高 适 生 区 的 总 面积 ;20 和 7D) 分别 为 1 时 段 第 i 个 单位 机 格 质 心 的 经 度 和 纬度 ; x(D、 y(D 分 别 为 
1 时 段 高 适 生 区 质心 的 经 度 和 纬度 。 


D= (x(r71)-3()) OCD- " 
RE - 
式 中 :DD 为 时段 到 1+1 时 段 高 适 生 区 的 位 移 距离 ;9 为 ?时 段 到 1+1 时 段 高 适 生 区 的 位 移 方向 , 0<9 <90° 


表示 位 移 方向 为 东北 , 909« 0 <180° 表 示 位 移 方向 为 西北 , 1809« Ó <270° 表 示 位 移 方向 为 西南 ， 
270°< 0 <360° 表 示 位 移 方向 为 东南 。 


2 结果 与 分 析 


2.1 影响 中 华 猕 猴 桃 分 布 的 关键 环境 变量 的 筛选 及 模拟 评价 
2.1.1 潜在 环境 变量 的 筛选 
本 研究 选择 的 初始 环境 变量 包括 预测 物种 分 布 通用 的 19 个 生物 气 RAE RURI 昌平 均 气 候 数据 ， 


Hj] 94-351. ix E 
程 中 会 引入 元 余 信 
Worthington 45:56 


研究 表 


ri ARN 避免 的 自 相 关 及 多 重 线性 重复 等 问题 ， 这些 高 度 相关 的 变量 在 模型 预测 过 
言 息 ， 影 响 预 测 结果 。 因 此 ， 应 首先 对 环境 变量 进行 相关 分 析 和 有 效 得 选 。 本 研究 中 , 参考 
方法 ， 筛 选 潜在 环境 变量 ， 有 具体 过 程 为 : 首先 全 部 环境 变量 构建 初始 模型 ， 选择 MaxEnt 软 


件 中 的 刀 切 法 (ackknife) 检 验 来 测定 环境 变量 对 模型 预测 的 贡献 大 小 ， 以 此 剔除 对 MaxEnt 模型 预测 结果 页 
献 较 小 的 环境 变量 。 其 次 , 利用 SPSS 软件 对 贡献 较 大 的 关键 限制 因子 进行 spearman 相关 分 析 ,， 若 相关 系 


"r-0.8, T 


共 保 留 了 22 个 环境 变量 (代码 和 计算 单位 见 表 2) 在 此 基础 上 重 


对 比 初始 模型 中 二 者 的 贡献 率 ,， 剔除 贡献 率 较 小 的 变量 ， 提 高 模型 模拟 的 精度 。 经 过 上 述 过 程 ， 


En 


E rp AE ROB E C rn E 746 BO e AA 78, 


并 对 模拟 结果 进行 准确 性 评价 。 
表 2 影响 中 华 猕猴 桃 分 布 的 22 个 环境 变量 及 其 代码 和 计量 单位 
Table 2 Codes and units of environmental variables used for simulation of potential distribution of Actinidia chinensis 
代码 变量 名 称 单位 
Code Environmental variable Unit 
Bio2 平均 日 较 差 Mean diurnal temperature range C 
Bio5 最 暖 月 最 高 温度 Max temperature of the warmest month C 
Bio6 最 冷 月 最 低温 度 Min temperature of the coldest month € 
Bio7 年 均 温 变化 范围 Annual temperature range ‘© 
Bio9 最 干 季 度 平均 温度 Mean temperature of the driest quarter jo 
Bioll 最 冷 季度 平均 温度 Mean temperature of the coldest quarter C 
Biol2 年 降水 量 Annual precipitation mm 
Biol4 最 干 月 降水 量 Precipitation of the driest month mm 
Precs, 9, 12 5. 90 和 12 月 平均 雨量 Precipitation in May, September, December mm 
2H. 4H. 9 H. 10 H. 11 月 和 12 月 最 高 温度 i 
Tmax2, 4, 9, 10, 11, 12 C 
Maximum temperature in February, April, October, November, December, and September 
Tmin3, 4, 10, 11 3H. 4H. 10 月 和 11 月 最 低温 度 Minimum temperature in March, April, October, and November T 
Tmean5 5 月 平均 温度 Mean temperature in May © 
2.1.2 ROC 曲线 和 AUC 值 对 模型 准确 性 评价 
图 1 是 当前 时 段 气候 条 件 下 的 ROC 曲线 图 , 图 中 表明 基于 主导 环境 变量 构建 的 中 华 猕猴 桃 在 中 国 的 地 


响 。 


0s. 


灵敏 度 (1- 遗 漏 率 ) 


Sensitive(1-omission ra 


AUC 平 均值 Mean(AUC=0.960) 日 
标准 误差 Mean+/-one stddev 
随机 预测 Random prediction — m 


00 0.1 02 03 04 05 06 07 08 09 10 
1- 特 异 率 (小 数 预测 面积 ) 
1-Specificity(fractional predicted area) 


图 1 当前 气候 情景 中 华 猕 猴 桃 MaxEnt 模型 的 ROC 曲线 


Fig.1 ROC curve of MaxEnt model for Actinidia chinensis under current scenario 


理 分 布 模型 的 AUC (873 0.960, 依据 表 1 的 AUC 值 评价 标准 ,本 次 构建 模型 的 预测 准确 性 达到 “ 极 好 ”的 标 
准 ; 图 2 为 未 来 不 同 气候 情境 下 MaxEnt 模拟 结果 的 ROC 曲线 图 , 由 图 可 知 所 有 预测 结果 的 AUC 
到 “ 极 好 ”的 标准 。 上 述 结果 说 明 模 型 预测 可 信和 度 较 高 ,可 用 于 分 析 气 候 变 化 对 中 华 猕 猴 桃 在 中 国 分 布 的 影 
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2030s, RCP2.6, AUC=0.958 , 2030s, RCP4.5, AUC=0.949 2030s, RCP8.5, AUC-0.952 


0 0.2 04 9.6 08 10 0 0.2 94 9.6 0.8 10 L] 92 9.4 9.6 98 10 


2050s, RCP2.6, AUC-0.964 2050s, RCP4.5, AUC-0.957 2050s, RCP8.5, AUC-0.955 


0 92 94 9.6 0.8 10 0 02 0.4 0.6 0.8 1.0 0 92 0.4 9.6 0.8 1.0 


2070s, RCP2.6, AUC-0.960 , 2070s, RCP4.5, AUC-0.956 , 2070s, RCP8.5, AUC-0.963 


0 02 94 9.6 0.8 1.0 0 02 0.4 0.6 0.8 1.0 0 92 0.4 9.6 0.8 1.0 


2080s, RCP2.6, AUC-0.954 , 2080s, RCP4.5, AUC-0.957 2080s, RCP8.5, AUC-0.964 


d 特异 率 (小 数 预 测 面积 ) s | 
1-Specificity(fractional predicted area) 
图 2 3 种 气候 变化 情景 下 中 华 猕猴 桃 MaxEnt 模型 的 ROC 曲线 


Fig.2 ROC curves of MaxEnt model for Actinidia chinensis under 3 climate change scenarios 


22 中华 猕 猴 桃 在 中 国 的 潜在 分 布 

由 预测 结果 (图 3) 和 统计 分 析 ( 表 4) 可 知 ， 当 前 气候 条 件 下 ,中华 猕猴 桃 在 中 国 的 适 生 区 在 25°~36?°N,， 
101°~122°E， 由 东 向 西 呈 不 连续 分 布 ， 总 面积 达 2.58x105km?， 占 中 国 国 土 面积 的 26.92%。 其 中 高 适 生 区 主 
要 位 于 四 川 、 陕西 、 重 庆 、 湖北 、 贵 州 、 浙 江 、 湖北 、 安 徽 、 河 南 、 江苏 和 甘肃 等 省 份 ， 面积 达 1.01x105km2， 
占 总 适 生 面积 的 38.94%; 中 适 生 区 沿 高 适 生 区 周围 分 布 ， 主 要 在 中 国 中 东部 地 区 , 包括 河南 、 湖 北 、 安 徽 、 


江苏 和 山东 等 地 ， 面 积 为 6.79x105km2， 占 总 适 生 面积 的 26.26%; 低 适 生 区 包括 云南 、 江 西 、 河 北 和 辽宁 等 
省 份 , 面积 为 8.99x105km2， 占 总 适 生 面积 的 34.896 
本 研究 选取 的 未 来 时 段 为 2030s、2050s、2070s 和 2080s, 3 种 气候 模式 为 RCP2.6、RCP4.5 和 RCP8.5, 


H 


JEEE: 


IK 4 可知， 上述 情 景 下 中 华 猕猴 桃 在 研究 区 的 高 、 中 适 生 区 仍 主要 集中 在 四 川 、 陕 西 、 重 庆 、 湖 


北 等 省 区 ， 只 是 在 不 同时 段 、 不 同 气候 情景 下 分 布 范围 及 面积 有 所 不 同 。 
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适 生 区 划 Suitability 
C] MEEK Unsuitable area 低 适 生 区 Lowly suitable area 


E 中 适 生 区 Moderately suitable area WEEE 高 适 生 区 Highly suitable area 


图 3 基于 MaxEnt 模型 预测 的 中 华 狗 猴 桃 在 中 国 的 适 生 分 析 图 


Fig.3 Potential distribution of Actinidia chinensis in China based on MaxEnt model 
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Eu 中 适 生 区 Moderately suitable area. BB =£ [X Highly suitable area 
图 4 未 来 3 种 气候 变化 背景 下 中 华 猕 猴 桃 潜在 适 生 区 的 预测 结果 


Fig.4 Predicted future suitable areas for Actinidia chinensis under 3 climate change scenarios 


表 4 _ 中华 猕 猴 桃 在 当前 情景 (1950 一 2000 年 ) 及 未 来 气候 条 件 下 的 适 生 区 面积 预测 


Table 4 Predicted suitable areas for Actinidia chinensis under current and future climatic conditions 


低 适 生 区 中 适 生 区 高 适 生 区 
. 气候 情景 Lowly suitable area Moderately suitable area Highly suitable area 
年 代 2 MA e hE EL y 2 MA e pE EL N y .x e hE E A 
i dimas yoa 丰 p Pd mage CU EUN 预测 面积 占 当 
scenarios Predicted area 从 Predicted ara 7 积 比 例 Predicted 积 比例 
(X10* km?) Proportion of current (X10 km?) Proportion of current area Proportion of current 
predicted area (96) predicted area (%) (x10^km?) predicted area (96) 
当前 (1950 一 2000 年 
i MERE 一 89.91 一 67.86 一 100.62 一 
21 世纪 30 年 代 RCP2.6 91.77 102.07 70.39 103.37 108.84 108.18 
2030s RCP4.5 01.27 12.64 86.81 127.92 86.89 86.36 
RCP8.5 122.19 135.90 62.01 91.38 106.94 106.29 
21 世纪 50 年 代 RCP2.6 02.34 13.83 66.93 98.62 108.43 107.77 
2050s RCP4.5 93.77 104.29 59.29 87.37 116.97 116.25 
RCP8.5 12.79 25.44 65.02 95.82 110.06 109.39 
21 世纪 70 年 代 RCP2.6 45.85 一 75.87 111.81 113.10 112.41 
2070s RCP4.5 11.01 23.55 79.71 117.15 108.49 107.83 
RCP8.5 20.82 34.38 58.53 86.24 122.72 121.96 
21 世纪 80 年 代 RCP2.6 89.82 99.91 81.41 119.96 108.28 107.62 
2080s RCP4.5 80.19 89.19 68.21 100.52 148.29 148.37 
RCP8.5 100.49 111.77 80.21 118.19 109.38 108.70 


2.5 ”中 华 猕 猴 桃 在 中 国 高 适 生 区 面积 变化 及 中 心 点 位 移 

当前 及 RCP2.6 情景 下 ， 高 适 生 区 面积 呈 波 动 增加 的 变化 趋势 。 总 体 上 由 当前 的 1.01x 106 km? 增加 至 
2080s 的 1.08x106km2( 表 4). 当前 至 2080s， 高 适 生 区 中 心 点 由 重庆 云 阳 ( 当 前 ) 经 重庆 巫山 (2030s)、 重 庆 奉 
节 (2050s 和 2070s) 位 移 至 湖 北 恩施 (2080s)( 图 5), 总 体 上 疝 东南 方 同 移动 82.39 kmC 5). 
当前 及 RCP4.5 情景 下 ， 高 适 生 区 面积 呈 先 减少 后 增加 的 变化 趋势 。 总 体 上 由 当前 的 1.01x 105 km? 增加 
Æ 2080s 的 1 48x106 kE 4)。 当 前 至 2080s， 高 适 生 区 中 心 点 由 重庆 云 阳 (当前 和 2030s) 经 湖北 恩施 
(2050s)、 重 庆 素 节 (2070s) 位 移 至 重庆 巫 溪 (2080s)( 图 5)， 总 体 上 向 西北 方向 移动 32.44 km( 表 5) 。 

当前 及 RCP8.5 情景 下 ， 高 适 生 区 面积 呈 先 增加 后 减少 变化 趋势 .总体 上 由 当前 的 1.01106 km? 增加 至 
2080s 的 1.09x106km2( 表 4)。 当 前 至 2080s， 高 适 生 区 中 心 点 由 重庆 云 阳 (当前 ) 经 重庆 巫山 2030s)、 重 庆 巫 
溪 (2050s)、 四 川 宣 汉 (2070s) 位 移 至 重庆 巫 溪 (2080s)( 图 5), 总 体 上 向 东北 方向 移动 40.82 kmC 5) 。 

R5 3 种 气候 变化 情景 下 中 华 猕 猴 桃 高 适 生 区 中 心 点 位 移 距 离 和 方向 
Table 5 Shift distance and direction of mean center of highly suitable area of Actinidia chinensis under 3 climate change scenarios 


RCP2.6 气候 情景 SRES-RCP2.6 RCP4.5 气候 情景 SRES-RCP4.5 RCP8.5 气候 情景 SRES-RCPS.5 


Eh 


时 期 位 移 S 自 度 位 移 zn 度 位 移 e E 
: : 方向 . 方向 . 方向 
Period Displacement Angle Displacement m Angle Displacement m» Angle 
Direction Direction Direction 
(kn) C) (km) (C) (km) (°) 
当前 至 2030s 东南 东南 东北 
108.84 279.2 9.68 290.23 106.56 66.37 
From current to 2030s Southeast Southeast Northeast 
2030s 至 2050s 西北 东南 西南 
58.65 94.87 47.70 344.09 73.01 265.52 
From 2030s to 2050s Northwest Southeast Southeast 
2050s 至 2070s 东北 东北 西北 
60.54 31.19 33.70 32.51 88.50 97.14 
From 2050s to 2070s Northeast Northeast Northwest 
2070s 至 2080s 东南 西北 东南 
62.73 357.39 74.00 132.39 92.80 281.38 
From 2070s to 2080s Southeast Northwest Southeast 
当前 至 2080s 东南 西北 东北 
82.39 316.83 32.44 153.71 40.82 43.31 


From current to 2080s Southeast Northwest Northeast 
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"- € 当前 情景 下 中 心 点 位 置 Location of mean center under current scenario 
= 
~ 全 一 >-RCP2.6 情 景 下 中 心 点 位 置 及 位 移 Location and displacement of mean centers under scenario RCP2.6 
© B — RCP4.5 情 景 下 中 心 点 位 置 及 位 移 Location and displacement of mean centers under scenario RCP4.5 
C3 
n e 十 RCP8.5 情 景 下 中 心 点 位 置 及 位 移 Location and displacement of mean centers under scenario RCP8.5 
o "T" 
A 5 3 种 气候 变化 情景 下 中 华 猕 猴 桃 高 适 生 区 中 心 点 迁移 轨迹 
CN Fig.5 Variation of mean centers of highly suitable areas under different 3 climate change scenarios 
FK 
rv 3 讨论 与 结论 
EN MaxEnt 软件 基于 最 大 焙 原理 ， 以 物种 分 布 变 量 和 环境 变量 为 基础 ,利用 数学 模型 统计 分 析 炉 最 大 时 物 


种 的 分 布 状态 ,众多 研究 表明 B71， 在 分 布 数据 较 少 时 MaxEnt 模型 的 预测 结果 比 同类 预测 模型 更 精确 ， 因 此 
本 研究 选取 MaxEnt 软 件 来 预测 并 分 析 中 华 猕猴 桃 的 适 生 区 及 适 生 等 级 。 目 前 普遍 应 用 ROC 曲线 法 ( 即 AUC 
ge 法 ) 对 模型 精度 进行 评价 B8, AUC 的 取 值 范围 在 0.5-1, 值 越 接 近 于 1 模型 精度 越 高 B9]。 本 研究 结果 表明 ,， 当 
c BU JR oe Up F PEIER ^E DC TL AUC 值 均 大 于 0.9， 预 测 准确 性 达 “ 极 好 ”标准 , 说明 此 次 
im T5 78 TL EIS] 8 7-8 E np EB Sz Ba 28 30.6 BEBE, RI AH F P EAD E F EL HER 738 E; UIROS RR 
FI 的 研究 。 另 外 , 本 研究 利用 地 理 信息 系统 软件 ArcGIS 对 MaxEnt 输出 的 栅 格 文件 进行 后 期 处 理 , 使 目标 物 
种 的 分 布 数据 与 环境 变量 数据 在 栅 格 单元 上 相对 应 ， 有 效 降低 了 系统 误差 ,进一步 提高 了 数据 的 准确 性 。 
本 研究 利用 ArcGIS 软件 统计 了 未 来 3 种 气候 变化 情景 下 中 华 猕 猴 桃 适 生 区 面积 变化 情况 ,结果 显示 ， 
至 2080s, 高 、 中 适 生 区 面积 在 RCP2.6、RCP4.5 和 RCP8.5 情景 下 均 有 所 增加 ， 只 是 增加 的 幅度 不 同 。 这 
说 明 在 不 同 排放 情景 下 气候 变化 对 中 华 猕猴 桃 潜在 分 布 的 影响 有 一 定 的 不 确定 性 。 环 境 变 量 的 选择 和 作物 
分 布点 数据 的 多 少 对 MaxEnt 模型 的 模拟 准确 性 有 较 大 影响 。 张 杰 等 9 利用 MaxEnt 分 析 了 中 等 温室 气体 
排放 情景 下 (RCP4.5) 中 华 狂 猴 桃 在 中 国 潜在 分 布 区 的 变化 趋势 ,结果 表明 ， 至 2050 年 ,中华 猕猴 桃 适 生 区 
面积 较 当 前 情景 略 有 下 降 , 此 结论 与 本 文 研究 结果 并 不 完全 吻合 ,可 能 是 选择 了 不 同 的 变量 和 不 同 的 分 布 
点 数据 所 致 。 
气候 是 决定 地 球 上 物种 分 布 的 最 主要 因素 ,而 物种 分 布 格局 的 变化 则 是 对 气候 变化 最 明确 和 直接 的 反 
应 。 气 候 变 暧 可 能 实质 性 地 改变 陆地 生态 系统 的 结构 和 功能 ,使 生物 栖息 地 范围 与 分 布 区 发 生变 化 。 本 文 
参考 Yue 等 B31 的 计算 方法 ， 以 面积 作为 权重 , 计算 了 不 同等 级 适 生 区 中 心 点 位 置 ， 着 重 分 析 了 适 生 区 中 心 
点 随时 间 变 化 的 位 移 趋势 ， 展示 了 中 华 猕 猴 桃 对 气候 变化 的 响应 过 程 ， 结 果 显 示 ， 到 2080s， m PEEK 
的 中 心 点 均 将 发 生 不 同 程度 的 改变 。 其 中 高 适 生 区 中 心 点 在 3 种 情景 下 均 向 偏 北方 向 移动 ， 原因 可 能 与 研 
究 区 的 年 平均 气温 与 年 降水 量 的 增加 有 关 ( 表 6); 从 结果 中 还 可 看 出 , 不 同 概率 等 级 的 中 心 点 位 没有 明显 的 
规律 性 ， 其 原因 可 能 在 于 当前 和 未 来 气候 变量 的 不 连贯 ,当前 气候 数据 的 起 止 时 间 为 1950—2000 年 ， 而 未 
来 气候 数据 的 起 止 时 间 为 2030—2080 年 , 2000 一 2030 年 的 数据 缺失 可 能 造成 中 心 点 移动 方向 的 往复 。 


表 6 气候 变化 情景 下 中 国 年 平均 气温 与 年 降水 量 的 变化 


Table6 Changes in annual mean temperature and annual precipitation under different climate change scenarios in China 


气候 变化 情景 年 份 年 平均 气温 范围 年 平均 气温 年 降水 量 范围 年 降水 量 

Climate change Feii Range of annual mean Annual mean Range of annual Annual precipitation 
scenario temperature (C) temperature (C) precipitation (mm) (mm) 

当前 Current 1950—2000 —13.8-25.8 6.15 12-4 199 621.59 

RCP2.6 气候 情景 2030s -18.2-27.1 6.68 16-5 062 824.16 

SRES-RCP2.6 2050s -17.7-21.3 6.94 18-5 149 857.32 

2070s -17.6-274 6.89 22-5 155 849.40 

2080s -17.9-27.3 6.79 22-5 101 845.79 

RCP4.5 气候 情景 2030s —18.0-27.2 6.73 16-5 086 822.52 

SRES-RCP4.5 2050s -17.1-27.6 7.28 19-5 201 849.65 

2070s —16.6-27.8 7.74 22~5 155 858.04 

2080s —16.4~27.9 7.69 21~5 143 859.39 

RCP8.5 气候 情景 2030s -17.5-274 11.92 21-5 054 819.23 

SRES-RCP8.5 2050s -16.4-27.9 11.65 21-5 149 858.42 

2070s -14.4-29.0 11.49 27-5 261 2024.12 

2080s -13.4-29.4 11.99 29-5 286 2 075.86 
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生态 位 模型 假设 物种 的 生态 位 需求 是 保守 的 , 样本 大 小 、 空 间 尺 度 和 环境 变量 等 因素 会 影响 此 类 模型 
的 预测 能 力 和 稳定 性 中 。 本 研究 中 发 现 以 下 问题 : 1) 文 中 选取 的 19 个 环境 变量 是 以 温度 和 雨量 为 基础 , 根 
据 不 同 需求 计算 演 生 而 来 , 因此 这 些 变量 之 间 存 在 不 可 避免 的 自 相 关 及 多 重 线性 重复 等 问题 ,在 模型 预测 
过 程 中 会 引入 元 余 信 息 , 影响 预测 结果 。 因 此 本 文 首 先 运用 刀 切 法 对 环境 变量 进行 相关 分 析 和 有 效 筛选 ， 
计算 各 因子 对 物种 分 布 的 贡献 率 大 小 ,剔除 贡献 率 较 小 的 环境 变量 ， 基 于 主导 环境 变量 重建 模型 ， 提 高 了 
预测 结果 的 准确 性 .2) 本 研究 预测 了 未 来 情景 下 中 华 猕猴 桃 的 适 生 区 变化 ,使 用 RCP2.6、RCP4.5 和 RCP8.5 
3 种 气候 模式 数据 ,需要 注意 的 是 , 预测 是 以 当前 中 华 猕 猴 桃 的 发 生 点 作为 分 布 数据 变量 , 这 就 可 能 会 忽略 
未 来 实际 分 布点 的 影响 , 不 可 避免 地 造成 系统 误差 。3) 影 响 中 华 猕猴 桃 在 中 国 种植 分 布 的 环境 因子 不 仪 仅 
包括 气候 条 件 、 土 壤 类 型 、 植 被 类 型 和 地 形 因子 、 品 种 类 型 、 人 类 活动 等 非 生物 因子 ,社会 经 济 结构 、 生 
产 技术 水 平等 因素 同样 会 对 其 分 布 产 生 重 要 影响 中 41。 基 于 上 述 原因 ， 可 以 推测 , MaxEnt 模型 预测 的 生态 
位 比 中 华 猕 猴 桃 所 占据 的 实际 生态 位 要 宽 。 本 研究 由 于 数据 限制 , 仅 考虑 了 温度 和 降水 这 两 类 气候 变量 对 
适 生 分 布 的 影响 , 在 下 一 步 工 作 中 , 还 应 注重 考量 各 种 因素 相互 作用 的 可 信和 表达 ， 以 改善 模型 的 预测 效果 。 
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